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We previously observed that an intrinsic staking fault 
shrunk through a glide of a Shockley partial disloca- 
tion terminating its lower end in a hard-sphere crys- 
tal under gravity coherently grown in (001) by Monte 
Carlo simulations [Mori et al., Molec. Phys. 105, 1377 
(2007)]; it was an answer to a one-decade long stand- 
ing question why the stacking disorder in colloidal crys- 
tals reduced under gravity [Zhu et al., Nature 387, 883 
(1997)]. Here, we present an elastic energy calculation; 
in addition to the self-energy of the partial disloca- 
tion [Mori et al, Prog. Theor. Phys. Suppl. 178, 33 
(2009)] we calculate the cross-coupling term between 
elastic field due to gravity and that due to a Shock- 
ley partial dislocation. The cross term is a increasing 
function of the linear dimension R over which the elastic 
field expands, showing that a driving force arises for the 
partial dislocation moving toward the upper boundary 
of a grain. 

Keywords: elasticity, Shockley partial dislocation, grav- 
ity, hard-sphere crystal 

1 Introduction 

Zhu et al. 1 in 1997 presented results of colloidal 
crystallization, which implied that the stacking dis- 
order in the colloidal crystals reduces due to gravity. 
They made such conclusion by comparison between 
results of a Space Shuttle experiment and those 
on the ground that whereas under micro gravity 
the colloidal crystals are of random hexagonal-close 
pack (rhcp) structure, under normal gravity they 
exhibit a mixture of rhcp and the face-centered cu- 
bic (fee) crystal. In a previous experiment the struc- 
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ture under normal gravity was fee [2]. This trend 
was supported by Kegel and Dhont [3] . Their result 
was, however, slightly different; the structure was 
faulted twinned fee under normal gravity. In their 
experiments the hard sphere (HS) nature was real- 
ized by using poly(methyl methacrylate) (PMMA) 
microspheres suspended in a hydrocarbon medium 
[4|. The difference in the particle size and disper- 
sion medium may results in the variety of the fi- 
nal states under the normal gravity. We point out, 
however, that the variety is nature. It is natural 
to understand that the non-uniqueness of the final 
states under gravity results from metastable equi- 
libria. Processes of disappearing stacking disorder 
must be those between certain metastable states. 

We performed Monte Carlo (MC) simulations of 
HSs at various gravitational conditions to elucidate 
the mechanism of the defect disappearance due to 
gravity. The stacking sequence of the hexagonal 
planes does not affect the density. So, the mecha- 
nism has been a hot question since Zhu et al. [1] . It 
was found in this simulation that a defective crys- 
tal was formed above a less-defective crystalline re- 
gion at the bottom of the system [5]. The lateral 
cross sections were square. Therefore, in case that 
the lateral system size was small, (001) planes were 
stacked along vertical direction. On the other hand, 
in case that the lateral system size was not so small, 
until the strength of gravity was not so large that 
the packing of particles at the bottom was not dom- 
inated by the smallness of the system (111) planes 
were stacked. One can understand those phenom- 
ena as results from competition between contribu- 
tions of the interfacial free energy against bottom 
wall and stress from the lateral boundary. We wish 
to point out that stress equivalent in those Simula- 
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tions can be given by use of a patterned substrate; 
the use of a patterned substrate was proposed by 
van Blaaderen et al. [6]. This method is called a 
colloidal epitaxy. We have taken a detailed look 
into processes in the (001) growth, in which a por- 
tion of the defective region is transformed into the 
less defective state 0. In the (001) growth {111} 
planes and, thus, stacking faults run along a oblique 
direction, not horizontal. 

From illustrations {e.g. Fig. 1) we have an intu- 
ition that the strain due to a Shockley partial dis- 
location at the lower end of an stacking fault coop- 
erates with gravitational effect toward defect disap- 
pearance. To give support to this intuition we have 
performed elastic theoretical calculations. A part of 
them has already been presented [8] ; the self-energy 
of the partial dislocation as well as the buoyancy 
due to particle deficiency of the core structure yields 
a driving force for the partial dislocation moving 
upward. Interfacial energy due to stacking disorder 
has been added as well as the dislocation core en- 
ergy. The former is a form of the diving force for 
shrinking of a staking fault running upward. Sec- 
tion [3] repeats some details of the elastic theoretical 
calculation presented already [8]. 

Strain energies due to a patterned substrate were 
evaluated for wetting by a crystalline layer [9l [10] . 
Also strain energy was estimated for misfit disloca- 
tion in colloidal epitaxy [lT| . The critical thickness 
he for misfit dislocation was estimated and for col- 
loidal crystals thicker than he a position where the 
stacking fault starts upward was given. Compar- 
ison between strain energies due to the substrate 
and due to lateral boundary would answer whether 
in phenomena in our MC simulation under small 
system size corresponds to those in colloidal epi- 
taxy. Strain in our MC simulation is geometrically 
equivalent to those which appear due to patterned 
substrate [12]. The interplay between gravitational 
and substrate effects is left for a future research. 
We will concentrate in this paper on the stability of 
an intrinsic stacking fault terminated by a Shockley 
partial dislocation. 

2 System 

An intrinsic stacking fault disappearing is schemat- 
ically shown in an illustration (Fig. 1). It is seen 
that a Shockley partial dislocation of the Burgers 
vector = (1/6)[211] (= -ai/3 + 02/6 + 03/6 
with tti, a2, and 03 being the lattice vectors) termi- 
nates the stacking fault. Translation of (111) plane 
by shifts the positions of particles in the (111) 



single Shockley partial viewed from [110] 




I [110] 

Figure 1: Illustration (not a snapshot of a simu- 
lation) of an intrinsic stacking fault which is disap- 
pearing. Particles in the distorted crystal are shown 
by the dots and the regular lattice positions by open 
circles. A dotted line indicates the stacking fault. In 
this illustration particles outside the portion right 
of this line are fixed at the regular positions. The 
intrinsic stacking fault is terminated by a Shockley 
partial dislocation at the center of this illustration. 

plane from A, B, and C to B, C, and A, respectively 
(see, Fig. 2.) Here, A, B, and C refer to three pos- 
sible positions in projection onto the (111) plane. 
Pictures similar to Fig. 1 obtained by confocal mi- 
croscopy have been presented by Schall et al. [11] 
and Alsayed et al. [T3]. Phenomena observed by 
Schall et al. are, however, different from what we 
focus in this paper. They have focused on misfit dis- 
locations and have observed increasing dislocations. 
On the other hand, Alsayed et al. observed premelt- 
ing near varieties of defects. Though the purposes 
are different, the Shockley partial dislocations evi- 
dently exist in colloidal crystals. The stacking fault 
runs upward and the linear dimension R over which 
the elastic field expands is regarded as the distance 
from the Shockley partial dislocation to the upper 
boundary of the crystalline grain. We will not con- 
sider the interaction between the bottom substrate 
and the Shockley partial dislocation. 

We note here that the density decreases with the 
altitude z, according to the mechanical equilibrium 
equation 

dP 

— = -mgp{z), (2.1) 
oz 

where P is the pressure, m the mass of a particle, g 
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Figure 2: The Burgers vector b' = (1/6)[211] (ar- 
row) and (111) plane (painted). The arrow connects 
a lattice position at, say, A and an adjacent one at, 
say, B. 

the acceleration due to gravity, and p{z) the particle 
number density at z in coarse scale. The other point 
characteristic of these simulations is that the coher- 
ent growth |12j occurred as already mentioned. We 
find that the lattice constant c is linearly increases 
with the altitude z whereas the lattice constants 
a — b are constant. We can write 



c(z) = c(0)[l + C^] 



(2.2) 



with ^ being the coefficient of order 10~^ [12]. In 
section d] we integrate fa^ ~ dua^/dxp to obtain 
the stress a[^p due to gravity with use of Eqs. (|2.ip 

and JMI), where fi^^ = dP/dxa- Here, the Greek 
alphabets, a, /?, have represented the Cartesian co- 
ordinate X, y, and z. We will follow this convention 
hereafter. 



Self energy, buoyancy, 
stacking fault energy 



and 



We shall reproduce the calculation of the sum of 
strain, stacking fault, and gravitational energies 
for the configurations illustrated in Fig. 1 per unit 
length (cr) depth perpendicular to the paper, where 
a is the HS diameter. In Ref. [8j the elastic strain 
energy for an intrinsic stacking fault was calcu- 
lated by using a formula W{R) = (^6^/47r)[cos^ 
sin^0/(l — iy)]\n{aR/b), where fi, v, and 6 are the 
shear modulus, the Poisson ratio, and the angle be- 
tween the Burgers vector b and the sense direction 
of the dislocation line, respectively, and R is the lin- 
ear dimension over which the elastic field expands 
and ro = b/a defines the radius of dislocation core 



with a ^ 0.1 14J (Do not confuse the coefficient 
a with the subscript a). For the intrinsic stack- 
ing fault we put b = b^ {\b^\ — a/y/6 with a being 



the fee lattice constant) and 9 
obtain 



tt/6 for W{R) to 



pa 

mi 



1 



1 - V 



In 



(3.1) 



We have Usf = ^dsf for the stacking fault energy 
of the intrinsic stacking fault, where 7i is the inter- 
facial energy of the intrinsic stacking fault and Isf 
is the effective length of the stacking fault. Igf is 
the same order of i?; a geometrical factor should be 
multiplied to R. Pronk and Frenkel [I5j calculated 
the stacking fault energy ^sf for the interfacial free 
energy between fee and hep crystals. 7^ is regarded 
as the same order of 7^ j . A factor should be mul- 
tiplied to 7s/. Introducing a factor C we can write 
the stacking fault energy as Usf = ClsfR- 

Let us estimate the gravitational energy. Parti- 
cle deficiency due to one Shockley partial disloca- 
tion corresponds to one-third lattice line in Fig. 1. 
A simple estimation is Ug — {l/2y/2)mgpa^R/3, 
where l/2\/2 comes from the number of particles 
per unit length of the lattice line [110] 0. 

The sum of Usf and Ug yields a linear term 

Usf + Ug = {Cisf + mgpa /6V2)R. (3.2) 

Because jsfO'^ /ksT calculated [13] is of order 10"'', 
Usf is negligible small as compared to Ug (in MC 
simulations [3 [T] [8j [12] we focus on the phenom- 
ena at mga/kBT 1), where ksT is the temper- 
ature multiplied by Boltzmann's constant. More- 
over, pa^/ksT = 50 - 100 [16l dT] [E] ; hence, Uei 
dominates. 

Energy due to the dislocation core should be 
added. It depends on the core structure. In particu- 
lar, in the present case the entropy term contributes 
through the change in the free volume or vibrational 
mode due to the dislocation core. Hence, the accu- 
rate evaluation of core energies should be made with 
simulations. Let us borrow the argument for solid 
states; the core energy is proportional to pb'^. In 
most cases the metal core energy Ucore is about ten 
percents of Uei- Let us write Ucore — Ppb^ with 
l3 being a certain constant of order less than unity 
(Do not confuse this coefficient with the subscript 
/3). Putting b ~ for the intrinsic stacking fault 
we have 

Ucore ^ Ppa^ /&. (3.3) 



^This factor has been laclting in Ref. [§]. Also, g has been 
missing after m there. 
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The elastic constants for the HS crystal were cal- 
culated by Frenkel and Ladd [16] by molecular dy- 
namic simulations, by Runge and Chester [17] by 
a MC simulation, and by Laird |18] by a density 
functional theory. The shear modulus jj, for the 
HS crystal ranges between 50 and 100 as already 
mentioned (in the unit of ksT/a^), which depends 
on the particle number density p. This range of 
fj. corresponds to pa^ = 1.06 — 1.13 {a/a = 1.55 
- 1.52) around the disappearing stacking disorder 
in the MC simulations [3 HI [H]. Stacking fauh 
energy was calculated by MC simulation by Pronk 
and Frenkel [15]. They calculated interfacial free 
energy between fee and hep at pa^ = 1.10 to be 
(26±6) xlO~^kBT/a^, as mentioned already. We 
neglect difference among fcc-hcp, intrinsic (such as 
...CAB- ABC...), and others as done 15 . We take 
a distance between the Shockley partial dislocation 
to the fluid region in the simulation or a linear di- 
mension of the grain in experiment for the effective 
length of the intrinsic stacking fault Isf as well as 
R. The "total" energy of the intrinsic stacking fault 
is 



U 



1 

96n 



1 



Clsf 



mgpa 
6%/2 
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R. 



pa 
(3.4) 



We note here that though the materials parameters 
for the HS system have been used, the formulas, 
Eqs. (|3. 1113. 41) . Eq. ((3?5t and those derived there- 
from, are valid for any materials as long as the lat- 
tice constant varies linearly in vertical direction. 

It is instructive to write the true total elastic en- 
ergy in terms of the elastic fields excluding the one 
which is forcing the coherent growth. 



U, 



1 



f (edge) 



, (edge) (screw) ^ Jg) 



(screw) 



2 I "^"/^ 



(edge) (edge) 



dV 



(screw) (screw) it r 
2 / <fl dV 



Aa)Aedge) 



dV + \j o^f/^ldV (3.5) 



where and Eap are, respectively, the stress and 
strain with superscripts {edge), {screw), and {g) in- 
dicating the edge and screw components of the elas- 
tic field due to the Shockley partial dislocation and 
that due to gravity, respectively. In the above ex- 
pansion diagonal relation between (screw) and (g) 
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Figure 3: The coordination system x\X2Xz for cal- 
culation of the elastic field due to the Shockley par- 
tial dislocation. The hexagonal planes are shown by 
thick horizontal lines. The intrinsic stacking fault is 
hatched. 2:3 axis is taken to be perpendicular to the 
hexagonal planes ([111]), axis along the disloca- 
tion line ([110]), and axis parallel to the stacking 
fauh ([112]). 

as well as {edge) has been used, which shall be con- 
firmed directly by seeing the component of the elas- 
tic fields shown in the preceding section. To note 
the fact that the sum of first two term in Eq. p.5|) 
equals U^i puts the present issue in relief; only the 
term / '^'^ap^at^'^'^ d^ depends simultaneously on {g) 
and {edge). The self-energy term \ J crlfp^apdV of 
the elastic field due to gravity does not depend on 
the position of the Shockley partial dislocation. 



4 Calculation of cross term 

At first, we write the elastic fields due to gravity 
and the Shockley partial dislocation. Integrating 



dxR 



(4.1) 



with the gravitational force /^^^ equaling the gradi- 
ent of the pressure P [as shown by Eq. (|2.1I) the x 
and y components vanish] we have 

aif) = ailHO) - nigp{0)-^ ln(l -I- ^z). (4.2) 

Here, we have used the relation p{z) = p{Q)/[l+£^z] 
[insertion of Eq. into p{z) = 4/a6c(z)]. The 

stress components a'^p except for a/? — zz vanish. 
Accordingly, we have 



(g) (edge),y 
"cxl3^cxl3 "'^ 



ji'^eif^^UV. (4.3) 



To calculate e'ff^'^^ as done previously [S] we intro- 
duce coordination system X1X2X3 shown in Fig. 3. 
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X3 axis is taken to be perpendicular to the hexag- 
onal plane, X2 axis parallel to the dislocation line 
((l/-\/2)[IlO]), and xi axis [within (111) and per- 
pendicular to X2] so as to the coordination system 
to be a right-handed system. Thus, we have unit 
vectors pointing to the axes as 

x\ = ^[112], af2 = ^[110], X3 = ^[111]- (4.4) 

In the frame of X1X2X3 the non-vanishing compo- 
nent of the strain due to the Shockley partial dislo- 
cation with Burgers vector 5^ = b| + b^j^ is given as 

m 

X3[{2\ + 3n)xl + fixl] 



{screw 



By calculation e 

Now, let us transform 



vanishes. 

Xl 



x 

y 



The 



X2 I into 
X3 / \ z 

inverse transform of x^ — daiXi is given by the 
transpose of aai because aai is orthogonal. 
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\bi\ x3[i2\ + ^i)xl-fixl] 
{\ + 2^i){xl + x^,)^ ' 
xi{xl 



271 



K\ 



4) 
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u) {xl+xD^ 



(screw) 



{screw) 



X3 



{screw) {screw) 
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47r xf -\ 

Ibjl X, 
Att xl+ i 



(4.5) 
(4.6) 
(4.7) 
(4.8) 
(4.9) 



{edge) 



is rewritten in terms of a;,y, and z as 

1 



(edge) _ 



18V37r(a;2 + z2)2 \ X + 2n 
-V2{2X + ll/i)x^ + (6A - 3n)x'^z 
(4A + 13m)z^^ 
V2x^ + Gx'^z 



-9V2fixz^ 



1 - I' 



where 6j| and bj_ are, respectively, the parallel 
(screw) and normal (edge) components of to the 
dislocation line. The stress components are read- 
ily obtained |14] . As mentioned previously we con- 
firm the orthogonal relation, i.e., u(':'^9e) ^(s^crew) _ ^ 
Here and hereafter, subscripts i,j represent 1,2,3. 
The orthogonal relation between (g) and (screw) 
is confirmed by coordinate transformation between 
xyz and X1X2X3. This coordinate transformation is 
also necessary to calculate / cr'f^el^^^'^^ dV . 

The matrix of the coordinate transformation 
from X1X2X3 system to xyz system is given by the 
inner product of the axis vectors x~a ■ Xi . The result 
is as 



-5\/2a;z^ + 2z^) I 



(4.13) 



Let us calculate / cri^z* ei'^f^'^'' dV substituting by 
Eqs. (|4.2p and (14.13^ . Integration will be performed 
in the cylindrical coordinate (r, 0) with the cylin- 
drical axis being y axis, r spans [ro,-R] and (/) does 
[0,27r]. At first, J ail* {0)eif'^^Uv = is shown by 
simple calculation. So, we have 



rilhi'f^''>dV 




1 





V6 \ 
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1 1 


2 







V6 


V3 / 

(4.10) 



which is the mere rotation in xz plane. We calculate 



(edge 



mffp(0)i J Hl+^z)eifs^^dV. 

(4.14) 

The integrations including ln(l -I- ^z) are shown in 
Appendix. We perform the integration for r, and 
then expanding with respect to ^ and integrate for 
(p. In expansion ^i? is treated as a small quantity. 
In Ref. [H] ^ is of the order of 10~^; so, for the 
grain with linear dimension R of several hundreds 
of particles this assumption is justified. 

Substituting by Eq. (HT^ in Eq. ({TIM and then 
using Eqs. ()A.15llA.18P we obtain 



as follows. 
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(4.11) 
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Figure 4: The coefRcients of i?^ (top) and i?** (bot- 
tom) of the cross-couphng term [Eqs. (|4.16p and 
(|4.17p ] are plotted against pa^ using values of elas- 
tic constants calculated in Refs. [16l [171 [18] . Both 
coefficients are shown to be positive. 

Let us look at the mangitude of the coefficients of 
and R"^. In the square bracket in Eq. ()4.15p the 
constant term is 



6A + 42/i 12 
X + 2fj, 1- 1 

and the coefficient of (^{R"^ + r§) 

14A + 68^i 4 



IS 



12{\ + 2fj,) 3(1-1^)' 



(4.16) 



(4.17) 



Quantities of Eqs. (|4.16p and (|4.17l) are plotted 
against pa^ in Fig. [4] It is shown that both co- 
efficients are positive, indicating that the cross- 
coupling term also yields a driving force for the 
Shockley partial dislocation moving toward the up- 
per boundary of the grain. 

5 Concluding remarks 

We have considered the system including an intrin- 
sic stacking fault with a Shockley partial dislocation 
terminating its lower end under gravity. In addi- 
tion to the self energy of the elastic field due to the 
Shockley partial dislocation, which has been shown 
to yield a driving force for the Shockley partial dis- 
location moving toward the upper boundary of the 
grain [S] , the cross-coupling term between the elas- 
tic field due to gravity and that due to the Shock- 
ley partial dislocation has been calculated. Driving 



force for the Shockley partial dislocation moving to- 
ward the upper boundary of the grain is also shown 
to arise from the cross term. 

The coefficient of Ini? in Eq. p.ip is, at most, 
5. Taking into account the denominator, the coeffi- 
cient of R^ in Eq. ()4.15p is a little less than this coef- 
ficient. It shows that the cross-coupling term gives 
a contribution, which is not negligible. In contrast, 
the coefficient of i?* is negligibly small because of 
the prefactor of Of course, the magnitude of 
the coefficient ^ depends on the strength of gravity 
(on ma and l/k^T through mga/ksT in natural 
sedimentation, and directly on g in centrifugal sed- 
imentation [ini [101 [H] ) ■ Under an extremely high 
g condition, however, the linear relation between c 
and z, Eq. (j2.2p . may not hold. This regime is, nev- 
ertheless, interesting because more enhancement of 
gravitational effect is expected. Consideration on 
this regime is left as a future research. 

We wish to emphasize that the simulation of col- 
loidal epitaxy is underway. The epitaxial growth 
has been confirmed when the lateral system size 
was increased doubly and doubly. It is an indi- 
cation that the stress causing the epitaxial growth 
is realizable. Not only simulations but also elas- 
tic theoretical study on the effect of the patterned 
substrate under gravity is also a future theme. 



A Calculations of integrals 



ln(l + ^z] 



dxdz 
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Integration for r can be done by applying 
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g integrals. 
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/ / sin^ 
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(A.4) 



ln(l + sin (j))dr = 



1 



-(1 + ^rsint 



C sin ( 

X ln(l + ^rsin^) — r. (A. 5) 
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Right-hand sides (RHS) of Eqs. (IA.1IIA.4I) are cal- f^'^ 2 , ■ 3 , , , /" 2 , ■ 4 



culated to be ^ Jo 12 _ 



X [(1 + ^i?sin0)ln(l +^i?sin0) / cos0sin(?!)(l + asint/)) ln(l + asin0)d0 

-(l + ^rosin0)ln(l+erosin0)]d0, (A.6) 

x'^z\n(l + fz) , , 1 Z"^'' 9 , =a cos (j) sin^ (j)d(p + — / cos(^sin^ 

- - / cos Jo 2 Jo 

X [(1 + ^i?sin0) ln(l + ^i?sin0) — ^ / cos 4> sin^ 4)d(f> + — cos^sin^ 

- sin (f>) \n{l+^ro sin (f>)]d(j), (A.7) ° ^ ° 



xz2ln(l+ez) 1 Z-^- 

axaz = — cos ffl sm 



= 0, (A.13) 
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- sin (f)) \n{l+^ro sin (t>)]d(l), (A.8) 
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(a;2 + 22)2 ^ ^ -—J sin'^0d0+— y sm 

X [(l + ei?sin0)ln(l+Ci?sin(/.) "^23^ „4 5^ " 

-(l + erosin0)ln(l+erosin0)]d0. (A.9) " ° + Y T ^ 12 T " ' ' 

Here, J^^ cos'^ 4>d(j) = 0, /g^'^ cos2 </) sin = 0, Accordingly, up to the third order in ^ the fohowing 

/o^''cos0sin2(/)d(/) = 0, and J^"" sin^ (l)d(l) = have results are obtained, 
already been used. Let us expand (l + asin0) ln(l+ /" x"^ ln(l + ^z) 

asin(/)) with a representing ^i? and S^tq. / , 2 , — 2^2"'^^'^^ ^ (A. 15) 

J ' ^ ) 



/, • ,M /-, • ,\ • , " • 2 , a; zln(l + ^z) , , 

(1 + asni(/)) ln(l + asin0) = asm0 + — sm cp / — — — ^-^-^dxdz = 

I J (a:^ + z^)^ 



(a:2 + z2)2 

g-sin^0+gsin*0.... (A.10) +i!(^^l_IoK, (A.16) 



-,4 



(The exact form is (1 + X)ln(l + X) - X = f xz^ Injl +Jz) ^^^^^ ^ 

y^OO i-X)" + ^ N y (3,2.^22)2 ^ ^ 



^00 (^OO^liN y (2;2+22)2 

Z^n=l n(n+l) V 3 9 2 

Integrals of cos'^ (p/ sin ip, cos^ cos (p sin </>, sin^ (p f ^ ^'^(^ + ^'^) — "^^(-^ " ^o)'"" 



(2 2\2 

Eq. (|A.10[) are shown to vanish by simple calcu- 5^'^ (i?'* — Vq 

lations. Also, integrals of these terms multiplied by 96 
sin'^ (p in the third term are turned to vanish. And, 
for cos'^ (/)/ sin (/i and cos sine/) integrals after mul- 
tiplying by sin^ <p in the second term and sin"* <p in 
the fourth term also vanish. Thus, we have 

cos3 , . 

/ (1 -I- asm(/)) ln(l -I- asm 

Jo sm (p 

= a / cos'^ (pd(p + -;r cos"^ (p sin 
Jo 2 Jq 

q3 ^27r ^4 

cos'^ o!) sin2 cpdcp H — - / cos'^ sm" 



(A.18) 



6 7o ^ ^ 12 

0, (A.ll) 



^27r 

/ cos^ 0(1 -|- asint/!)) ln(l + asin 
Jo 

= a cos2 sin H / cos2 6 sin^ 
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